Evaluating the locality of intrinsic precession damping in transition metals 
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The Landau-Lifshitz-Gilbert damping parameter is typically assumed to be a local quantity, in- 
dependent of magnetic configuration. To test the validity of this assumption we calculate the 
precession damping rate of small amplitude non-uniform mode magnons in iron, cobalt, and nickel. 
At scattering rates expected near and above room temperature, little change in the damping rate is 
found as the magnon wavelength is decreased from infinity to a length shorter than features probed 
in recent experiments. This result indicates that non-local effects due to the presence of weakly 
non-uniform modes, expected in real devices, should not appreciably affect the dynamic response of 
the element at typical operating temperatures. Conversely, at scattering rates expected in very pure 
samples around cryogenic temperatures, non-local effects result in an order of magnitude decrease 
in damping rates for magnons with wavelengths commensurate with domain wall widths. While this 
low temperature result is likely of little practical importance, it provides an experimentally testable 
prediction of the non-local contribution of the spin-orbit torque-correlation model of precession 
damping. None of these results exhibit strong dependence on the magnon propagation direction. 
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Magnetization dynamics continues to be a techno- 
logically important, but incompletely understood topic. 
Historically, field induced magnetization dynamics have 
been described adequately by the phenomenological 
Landau-Lifshitz (LL) equation 

rii = — ItmIhi X H -f Aril x (m x H) , (1) 

or the mathematically equivalent Gilbert form 0, Q. 
Equation [T] accounts for the near equilibrium dynamics 
of systems in the absence of an electrical current, jm 
is the gyromagnetic ratio and A is the phenomenological 
damping parameter, which quantifies the decay of the 
excited system back to equilibrium. The LL equation 
is a rather simple approximation to very intricate dy- 
namic processes. The limitations of the approximations 
entering into the LL equation are likely to be tested by 
the next generation of magnetodynamic devices. While 
many generalizations for the LL equation are possible, we 
focus on investigating the importance of non-local contri- 
butions to damping. It is generally assumed in both ana- 
lyzing experimental results and in performing micromag- 
netic simulations that damping is a local phenomenon. 
While no clear experimental evidence exists to contradict 
this assumption, the possibility that the damping is non- 
local - that it depends, for example, on the local gradient 
of the magnetization - would have particular implications 
for experiments that quantify spin-current polarization 
for storage [5|] and logic ^ devices based on using 
this spin-current to move domain- walls, quantifying vor- 
tex and mode [3| dynamics in patterned samples, and 
the behavior of nano-contact oscillators 0, . 

While several viable mechanisms have been proposed 
to explain the dam pin g process in different systems 
[ni,[l2,[ll,[li,[ll,[lllI3, we restrict the scope of this pa- 
per to investigating the degree to which the assumption 



of local damping is violated for small amplitude dynam- 
ics within pure bulk transition metal systems where the 
dominant source of damping is the intrinsic spin-orbit in- 
teraction. For such systems, Kambersky's 14] spin-orbit 
torque-correlation model, which predicts a decay rate for 
the uniform precession mode of 
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has recently been demonstrated to account for the major- 
ity of damping [H, [l^. The matrix elements |r~„(k)p 
represent a scattering event in which a quantum of the 
uniform mode decays into a single quasi-particle electron- 
hole excitation. This annihilation of a magnon raises the 
angular momentum of the system, orienting the magneti- 
zation closer to equilibrium. The excited electron, which 
has wavevector k and band index m, and the hole, with 
wavevector k and band index n, carry off the energy and 
angular momentum of the magnon. This electron-hole 
pair is rapidly quenched through lattice scattering. The 
weighting function Wnm(M) measures the rate at which 
the scattering event occurs. The very short lifetime of 
the electron-hole pair quasiparticle (on the order of fs at 
room temperature) introduces significant energy broad- 
ening (several hundred meV). The weighting function, 
which is a generalization of the delta function appearing 
in a simple Fermi's golden rule expression, quantifies the 
energy overlap of the broadened electron and hole states 
with each other and with the Fermi level. 

Equation [51 which has been discussed extensively 
[13, EM [H, Eyl , considers only local contributions to the 
damping rate. Non-local contributions to damping may 
be studied through the decay of non-uniform spin- waves. 
Although recent efforts have approached the problem of 



2 



non-local contributions to the dissipation of non-coUinear 
excited states [U the simple step of generalizing 
Kambersky's theory to non-uniform mode magnons has 
not yet been taken. We fill this obvious gap, obtaining a 
damping rate of 
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for a magnon with wavevector q. We test the impor- 
tance of non-local effects by quantifying this expression 
for varying degrees of magnetic non-coUinearity (magnon 
wavevector magnitude). The numerical evaluation of 
Eq. [3] for the damping rate of finite wavelength magnons 
in transition metal systems, presented in Fig. [1] and the 
ensuing physical discussion form the primary contribu- 
tion of this paper. We find that the damping rate ex- 
pected in very pure samples at low temperature is rapidly 
reduced as the magnon wavevector |q| grows, but the 
damping rate anticipated outside of this ideal limit is 
barely affected. We provide a simple band structure ar- 
gument to explain these observations. The results are 
relevant to systems for which the non-coUinear excita- 
tion may be expanded in long wavelength spin- waves, 
provided the amplitude of these waves is small enough to 
neglect magnon-magnon scattering. 

Calculations for the single-mode damping constant 
(Eq. [3]) as a function of electron scattering rate are pre- 
sented in Fig. [1] for iron, cobalt, and nickel. The Gilbert 
damping parameter a — A/7M is also given. Damp- 
ing rates are given for magnons with wavevectors along 
the bulk equilibrium directions, which are (100) for Fe, 
(0001) for Co, and (111) for Ni. Qualitatively and quan- 
titatively similar results were obtained for other magnon 
wavevector directions for each metal. The magnons re- 
ported on in Fig.[T]constitute small deviations of the mag- 
netization transverse to the equilibrium direction with 
wavevector magnitudes between zero and 1 % of the Bril- 
louin zone edge. This wavevector range corresponds to 
magnon half- wavelengths between infinity and 100 lat- 
tice spacings, which is 28.7 nm for Fe, 40.7 nm for Co, 
and 35.2 nm for Ni. This range includes the wavelengths 
reported by Vlaminck and Bailleul in their recent mea- 
surement of spin polarization 4]. 

Results for the three metals are qualitatively similar. 
The most striking trend is a dramatic, order of magnitude 
decrease of the damping rate at the lowest scattering rate 
tested as the wavevector magnitude increases from zero 
to 1 % of the Brillouin zone edge. This observation holds 
in each metal for every magnon propagation direction 
investigated. For the higher scattering rates expected in 
devices at room temperature there is almost no change 
in the damping rate as the magnon wavevector increases 
from zero to 1 % of the Brillouin zone edge in any of the 
directions investigated for any of the metals. 

To understand the different dependences of the damp- 
ing rate on the magnon wavevector at low versus high 
scattering rates we first note that the damping rate 
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FIG. 1: Damping rates versus scattering rate. The preces- 
sion damping rates for magnons in iron, cobalt, and nickel 
are plotted versus electron scattering rate for several magnon 
wavevectors. A dramatic reduction in damping rate is ob- 
served at the lowest scattering rates. The Landau-Lifshitz A 
(Gilbert a) damping parameter is given on the left (right) 
axes. Electron scattering rate is given in eV on the top axis. 
Magnon wavevector magnitudes are given in units of the Bril- 
louin zone edge and directions are as indicated in the text. 



(Eqs. [2] & [3]) is a convolution of two factors: the torque 
matrix elements and the weighting function. The ma- 
trix elements do not change significantly as the magnon 
wavevector increases, however, the weighting function 
can change substantially. The weighting function 

Vl/"„„j(k,k q) « A„,k(eF)^m,k+q(eF) (4) 
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FIG. 2: Partial band structure of bcc iron. The horizontal 
black line indicates the Fermi level and the shaded region 
represents the degree of spectral broadening. The solid dot is 
a hypothetical initial electron state while the open circle is a 
potential final scattering state. (Initial and final state wave- 
vector separations are exaggerated for clarity of illustration.) 
The intraband magnon decay rate diminishes as the energy 
separation of the states exceeds the spectral broadening. 



contains a product of the initial and final state electron 
spectral functions 
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which are Lorentzians in energy space. The spectral func- 
tion for state \n, k), which has nominal band energy e„.k, 
is evaluated within a very narrow range of the Fermi level 
ep- The width of the spectral function h'^ is given by 
the electron scattering rate 7 = l/2r where r is the 
orbital lifetime. (The lifetimes of all orbital states are 
taken to be equal for these calculations and no specific 
scattering mechanism is implied.) The weighting func- 
tion restricts the electron-hole pair generated during the 
magnon decay to states close in energy to each other and 
near the Fermi level. For high scattering rates, the elec- 
tron spectral functions are significantly broadened and 
the weighting function incorporates states within an ap- 
preciable range (several hundred meV) of the Fermi level. 
For low scattering rates, the spectral functions are quite 
narrow (only a few meV) and both the electron and hole 
state must be very close to the Fermi level. 

The second consideration useful for understanding the 
results of Fig. [1] is that the sum in Eqs. [5] & [3] can be 
divided into intraband {n = m) and interband {n ^ m) 
terms. For the uniform mode, these two contributions 
correspond to different physical processes with the intra- 
band contribution dominating at low scattering rates and 
the interband terms dominating at high scattering rates 
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FIG. 3: Intraband and interband damping contributions in 
iron. The intraband and interband contributions to the damp- 
ing rate of magnons in the (100) direction in iron are plot- 
ted versus scattering rate for several magnitudes of magnon 
wavevector. Magnitudes are given in units of the Brillouin 
zone edge. 



For intraband scattering, the electron and hole occupy 
the same band and must have essentially the same energy 
(within hj). The energy difference between the electron 
and hole states may be approximated as f:n,k+ci — ^n,k ~ 
q • 9e„,k/9k. The generation of intraband electron-hole 
pairs responsible for intraband damping gets suppressed 
as q-i9e„_k/c?k becomes large compared to hj. Unless the 
bands are very flat at the Fermi level there will be few lo- 
cations on the Fermi surface that maintain the condition 
q ■ 9e„^k/c?k < h'f for low scattering rates as the magnon 
wavevector grows. (See Fig.[2])- Indeed, at low scattering 
rates when is only a few meV, Fig. [3] shows that the 
intraband contribution to damping decreases markedly 
with only modest increase of the magnon wavevector. 
Since the intraband contribution dominates the inter- 
band term in this limit the total damping rate also de- 
creases sharply as the magnon wavevector is increased 
for low scattering rates. For higher scattering rates, the 
electron spectral functions are sufficiently broadened that 
the overlap of intraband states does not decrease appre- 
ciably as the states are separated by finite wavevector 
(q • den,k/d'k < fi-y generally holds over the Fermi sur- 
face). Therefore, the intraband contribution is largely in- 
dependent of magnon wavevector at high scattering rates. 

The interband contribution to damping involves scat- 
tering between states in different bands, separated by the 
magnon wavevector q. Isolating the interband damping 
contribution reveals that these contributions are insensi- 
tive to the magnon wavevector at higher scattering rates 
where they form the dominant contribution to damp- 
ing (see Fig. [3]). To understand these observations we 
again compare the spectral broadening ?i7 to the quasi- 
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quasiparticle energy difference may be approximated as 



particle energy difference A^^ k — Cm.k+q — fn.k- The 
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A™J^ + 1 ■ ^'^^'w/^^- The interband energy spacings 
are effectively modulated by the product of the magnon 
wavevector and the slopes of the bands. At high scatter- 
ing rates, when the spectral broadening exceeds the ver- 
tical band spacings, this energy modulation is unimpor- 
tant and the damping rate is independent of the magnon 
wavevector. At low scattering rates, when the spec- 
tral broadening is less than many of the band spacings, 
this modulation can alter the interband energy spacings 
enough to allow or forbid generation of these electron- 
hole pairs. For Fe, Co, and Ni, this produces a modest 
increase in the interband damping rate at low scattering 
rates as the magnon wavevector increases. However, this 
effect is unimportant to the total damping rate, which 
remains dominated by the intraband terms at low scat- 
tering rates. 

Lastly, we describe the numerical methods employed 
in this study. Converged ground state electron densities 
were first obtained via the linear-augmented-plane-wave 
method. The Perdew-Wang functional for exchange- 
correlation within the local spin density approximation 
was implemented. Many details of the ground state den- 
sity convergence process are given in [23| . Densities were 
then expanded into Kohn-Sham orbitals using a scalar- 
relativistic spin-orbit interaction with the magnetiza- 
tion aligned along the experimentally determined mag- 
netocrystalline anisotropy easy axis. The Kohn-Sham 
energies were artificially broadened through the ad hoc 
introduction of an electron lifetime. Matrix elements of 
the torque operator T~ = [a~ ,Ti.so] were evaluated sim- 
ilarly to the spin-orbit matrix elements {a^ is the 
spin lowering operator and Hgo is the spin-orbit Hamil- 
tonian.) The product of the matrix elements and the 
weighting function were integrated over fc-space using the 
special points method with a fictitious smearing of the 
Fermi surface for numerical stability. Convergence was 



obtained by sampling the full Brillouin zone with 160'^ 
fc-points for Fe and Ni, and 160^ x 91 points for Co. 

In summary, we have investigated the importance of 
non-local damping effects by calculating the intrinsic 
spin-orbit contribution to precession damping in bulk 
transition metal ferromagnets for small amplitude spin- 
waves with finite wavelengths. Results of the calculations 
do not contradict the common-practice assumption that 
damping is a local phenomenon. For transition metals, 
at scattering rates corresponding to room temperature, 
we find that the single-mode damping rate is essentially 
independent of magnon wavevector for wavevectors be- 
tween zero and 1 % of the Brillouin zone edge. It is not 
until low temperatures in the most pure samples that 
non-local effects become significant. At these scatter- 
ing rates, damping rates decrease by as much as an or- 
der of magnitude as the magnon wavevector is increased. 
The insensitivity of damping rate to magnon wavevector 
at high scattering rates versus the strong sensitivity at 
low scattering rates can be explained in terms of band 
structure effects. Due to electron spectral broadening at 
high scattering rates the energy conservation constraint 
during magnon decay is effectively relaxed, making the 
damping rate independent of magnon wavevector. The 
minimal spectral broadening at low scattering rates - 
seen only in very pure and cold samples - greatly restricts 
the possible intraband scattering processes, lowering the 
damping rate. The prediction of reduced damping at low 
scattering rates and non-zero magnon wavevectors is of 
little practical importance, but could provide an accessi- 
ble test of the torque-correlation model. Specifically, this 
might be testable in ferromagnetic semiconductors such 
as (Ga,Mn)As for which many spin-wave resonances have 
been experimentally observed at low temperatures [25| . 
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